Phenomenology of Gravitational Aether 
as a solution to the Old Cosmological Constant Problem 
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One of the deepest and most long-standing mysteries in physics has been the huge discrepancy 
between the observed vacuum density and our expectations from theories of high energy physics, 
which has been dubbed the Old Cosmological Constant problem. One proposal to address this 
puzzle at the semi-classical level is to decouple quantum vacuum from space-time geometry via a 
modification of gravity that includes an incompressible fluid, known as Gravitational Aether. In this 
paper, we discuss classical predictions of this theory along with its compatibility with cosmological 
and experimental tests of gravity. We argue that deviations from General Relativity (GR) in this 
theory are sourced by pressure or vorticity. In particular, the theory predicts that the gravita- 
tional constant for radiation is 33% larger than that of non-relativistic matter, which is preferred 
by (most) cosmic microwave background (CMB), Ly-a forest, and 7 Li primordial abundance obser- 
vations, while being consistent with other cosmological tests at ~ 2a level. It is further shown that 
all Parametrized Post-Newtonian (PPN) parameters have the standard GR values aside from the 
anomalous coupling to pressure £4, which has not been directly measured. A more subtle prediction 
of this model (assuming irrotational aether) is that the (intrinsic) gravitomagnetic effect is 33% 
larger than GR prediction. This is consistent with current limits from LAGEOS and Gravity Probe 
B at ~ 2a level. 



I. INTRODUCTION 



The discovery of recent acceleration of cosmic expan- 
sion was one of the most surprising findings in mod- 
ern cosmology jj], ||. The standard cosmological model 
(also known as the concordance model) drives this ex- 
pansion with a cosmological constant (CC). While the 
CC is consistent with (nearly) all current cosmological 
observations, it requires an extreme fine-tuning of more 
than 60 orders of magnitude, known as the cosmological 
constant problem Q. More precisely, a covariant regu- 
larization of the vacuum state energy of a Quantum Field 
Theory (QFT), if it exists, acts just like the CC in linear 
order, but has a value many orders of magnitude larger 
than what is inferred from observations. 

If the QFT prediction of the cosmological constant is 
considered reasonable (and in lieu of an extreme fine- 
tuning), there is no choice but to abandon the idea 
that vacuum energy should gravitate. This, however, re- 
quires modifying Einstein's theory of gravity, in which all 
sources of energy gravitate. Attempts in this direction 
have been proposed in the context of massive gravity [|] , 
or braneworld models of extra dimensions such as cascad- 
ing gravity [0, g|, or supersymmetric large extra dimen- 
sions (e.g., Rj). However, efforts to find explicit cosmo- 



logical solutions that de-gravitate vacuum have proved 
difficult (e.g., g |). 

In , one of us proposed a novel approach to modified 
gravity in which the QFT vacuum quantum fluctuations 
(of linear order in the metric) are decoupled from grav- 
ity through the introduction of an incompressible perfect 
fluid called the Gravitational Aether. In this model, the 
right hand side of the Einstein field equation is modified 
as: 
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where G is the (only) constant of the theory and T' is 
the aether fluid which has pressure p' and four-velocity 
u'. Our metric signature is ( — ,+,+,+). Aether is con- 
strained by requiring the conservation of the energy- 
momentum tensor !),„, and the Bianchi identity: 



V M T' — -V T 



(3) 



which can be written in a similar form to the relativistic 
hydrodynamic equations: 
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where = u' ■ V, and Vi is the gradient normal to u' 
four-vector. Notice that Eqs. (|[^) can be combined to 
find a parabolic equation for the evolution of u', which 
generically has a well-defined initial value problem, at 
least for a finite time [04j . 

This modification of Einstein equations (|l|-||), if self- 
consistent and in agreement with other experimental 
bounds on gravity, could potentially constitute a solution 
to the old cosmological constant problem. We will show 
in this paper that none of these experimental bounds, as 
of yet, rule out this theory (at ~ 2er level) and that it is 
surprisingly similar to general relativity |45j] . 

Nevertheless, the new cosmological constant problem, 
i.e. the present-day acceleration of cosmic expansion is 
not addressed by the original gravitational aether pro- 
posal. In 11, [12], it is argued that quantum gravity 
effects in the presence of astrophysical black holes can 
naturally explain this phenomenon. In this proposal, the 
formation of black holes leads to a negative aether pres- 
sure, that is set by the horizon temperature of the black 
holes. However, in the present work we only focus on 
phenomenological implications of the classical gravita- 
tional aether scenario, and defer study of black hole-dark 
energy connection, which could be potentially very im- 
portant on cosmological scales at late times. Instead, we 
use a standard cosmological constant to model the late- 
time acceleration of cosmic expansion. Throughout the 
paper we set the speed of light c=l. 



II. COSMOLOGICAL CONSTRAINTS ON 
GRAVITATIONAL AETHER 



If the energy- momentum tensor of matter, T^ Vl can be 
approximated by a perfect fluid with constant equation 
of state p = wp and four- velocity m m , direct substitution 
into Eq. (g) shows that if: - •* - (1+t " )(3 " J ~ 1) 



Ua,P 



p, 



then the solutions to the gravitational aether theory are 
identical to those of General Relativity (GR) with a 
renormalized gravitational constant: 



G N ^G eS = (l + w)G 



N: 



(6) 



where Gn — 3G'/4. In other words, the gravitational 
coupling is not a constant anymore, and can change sig- 
nificantly for fluids with relativistic pressure. Not sur- 
prisingly, for vacuum equation of state w = — 1, G c s — 0, 
which implies that vacuum does not gravitate. 

In particular, in the case of homogeneous FLRW cos- 
mology where the perfect fluid approximation is valid, 
this theory predicts that the effective G that relates ge- 
ometry to the matter density p in Friedmann equation is 
different in the matter and radiation eras: 



Gn_ = G cS {w = 0) = 3 
G R ~ G cS (w = 1/3) 4" 



(7) 



This is the first cosmological prediction of this the- 
ory: radiation energy gravitates more strongly than non- 



relativistic matter. The expansion history in the radia- 
tion era depends on the product Gp ra d, and is constrained 
through different observational probes. The constraints 
are often described as the bound on the effective num- 
ber of neutrinos N u ^ e g, which quantifies the total ra- 
diation density p ra d- However, assuming only photons 
(that are constrained by CMB observation) and three 
neutrino species, with no more light particles left over 
from the very early universe, we can translate the con- 
straints to those on G e s by requiring G gp ia d(N v = 3) = 
GjvPrad(^ = 3 + AN V ). In particular, based on stan- 
dard thermal decoupling of neutrinos, Eq. (0) can be 
translated to AN U = 2.5, at least for a homogeneous 
universe Jill . Using this correspondence, we can now dis- 
cuss cosmological constraints on the gravitational aether 
scenario. 
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FIG. 1: Allowed regions with 2 a lines for D/H, Y p and 
'Li/H are shown. The upper and lower horizontal 
dashed lines indicate GR and gravitational aether 

predictions, respectively. The thickness of Y p means the 
uncertainty in measurements of neutron 

lifetime |T|, 14 . We can translate the vertical axis into 

AN U by using a relation G N /G R . ^ 1/(1 + 0.135AiY„). 



A. Big Bang Nucleosynthesis 

It has been known that the increase of the gravita- 
tional constant at around T = 0(1) MeV epoch induces 
earlier freezeout of the neutron to proton ratio because 
of a speed-up effect of the increased cosmic expansion. 
This raises the abundance of 4 He sensitively and deu- 
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terium (D) mildly, and can lower the abundance of 7 Be 
through 7 Be (n,p) 7 Li(p, a) 4 He (Note that the second p is 
thermal proton) . For a relatively large baryon to photon 
ratio r\ > 3 x 10~ 10 , the dominant mode to produce 7 Li 
is the electron capture of 7 Be at a later epoch through 
7 Be + e~ — > 7 Li + v e . Therefore, the decrease of 7 Be 
makes the fitting better because so far any observational 
7 Li abundances have been so low that they could not 
have agreed with theoretical prediction in Standard BBN 
(SBBN) at better than 3 a ||. 

In this study, we adopt the following observational light 
element abundances as primordial values: the mass frac- 
tion of 4 He, Y p = 0.2561 ± 0.0108 (stat) @, the deu- 
terium to hydrogen ratio, D/H=(2.80±0.20) x 10~ 5 @, 
and the 7 Li to hydrogen ratio Log 10 ( 7 Li/H) = —9.63 ± 
0.06 |lf| Wh ' Theoretical errors come from experimen- 
tal uncertainties in cross sections [fBl |l9l |2(i|] and neutron 
lifetime 1R 
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Comparing theoretical prediction with the observa- 
tional light element aubndances provides a constraint 
on Gn/Gr. Fig. I] shows the results of a comprehen- 
sive analysis for 'Tie, D, and 7 Li. We also plotted 
a band for baryon to photon ratio, r\ which was re- 
ported from CMB observations by WMAP 7- year, 77 = 
(6.225±0.170)xlO- 10 incaseofGj V /G J R = l [||. Then 
we can see that every light element agrees with the Grav- 
itational Aether theory within 2 a. It is notable that 7 Li 
in this theory fits the data better than that in SBBN. 
Performing \ 2 fitting for three elements with three de- 
gree of freedom, however, the model is allowed only at 
99.7% (3 a) in total. 

However, notice that the main discrepancy is with deu- 
terium abundance observed in quasar absorption lines, 
which suffer from an unexplained scatter. Moreover, deu- 
terium could be depleted by absorption onto dust grains 
that would make its primordial value closer to our pre- 
diction (see p2| for a discussion). 



B. Power Spectrum of Cosmological Fluctuations 

The Gravitational Aether theory can also be tested 
by considering the power spectrum of the CMB, just 
as a number of publications have recently investigated 
the apparent preference for extra relativistic degrees of 
freedom (see e.g., p3|-p5|). Using a modified version of 
Cmbeasy |^6|, |7|, we compute constraints on Gn/Gr 
from scalar perturbations in a scenario with three mass- 
less neutrino species (details are discussed in Appendix 
0). The 7-year CMB data from WMAP @ together 
with small-scale observations from the Atacama Cosmol- 
ogy Telescope (ACT) || yield G N /G R = 0.73tgj} 
at 95%-confidence. Just like any additional relativistic 
component can be compensated by a higher fraction of 
dark matter in order to keep the time of matter-radiation 
equality constant, there is a high amount of degeneracy 
between Gn/Gr and fl m h 2 and h (see Fig. |]). Recent 
data from the South Pole Telescope (SPT), which mea- 





Gn/Gr 


WMAP+ACT 


0.73±g;gi 


WMAP+ACT+SPT 


0.88«;i 7 


WMAP+ACT+Hubble+BAO+Sne 


0.89t°;i? 


WMAP+ACT+SPT+Hubble+BAO+Sne 




WMAP+ACT+Sne+Ly-a (free Y p ) 


0.68t° o f 5 


WMAP+ACT+SPT+Sne+Ly-a (free Y p ) 


0.90«J 7 



TABLE I: Summary of the constraints on Gn/Gr 
and the associated 95% confidence intervals for 
different combinations of observational data. 



sured the CMB power spectrum in the multipole range 
650 < I < 3000 significantly tightens the constraint and 
yields 0.881^ (for the combination of ACT and SPT 
data we have adopted the SPT treatment of foreground 
nuisance parameters). A similar effect can be seen when 
adding Baryonic Acoustic Oscillations (BAO) |2j| and 
constraints on the Hubble rate. Here we adopt the value 
of H = 73.8 ± 2.4 km" 1 Mpc" 1 @. Then, by break- 
ing the degeneracy between the matter content and h, 
the combination WMAP+ACT+BAO+Snc+Hubble re- 
sult in Gn/Gr, = 0.89^0^. The supernovae data of 
the Union catalog |3l| do not significantly contribute to 
this constraint. We note that both cases, i.e. adding ei- 
ther SPT data or adding the Hubble constraints to the 
basic WMAP+ACT set, move the gravitational Aether 
value of Gn/Gr — 0.75 to the border or just outside of 
the 95% confidence interval, while the General Relativity 
value of Gn/Gr, = 1.0 is well compatible with all com- 
binations of data. Consequently, the full combination 
of data, i.e. WMAP+ACT+SPT+Hubble+BAO+Sne, 
constrains Gn/Gr to 0.941^°. 

In contrast, observational constraints at lower red- 
shifts, in particular data of the Ly-a forest |32| prefer the 
Aether prediction. Furthermore, additional degeneracies 
with e.g. the Helium mass fraction Y p might shift the 
preferred values. Combining WMAP+ACT+Sne with 
the Ly-a forest constraints yields, Gn/Gr = 0.68jiQ'25 
at 95% level, with Yp as a free parameter. However, we 
should note that this result is more prone to systematic 
uncertainties due to the quasilinear nature of the Ly-a 
forest. Also, including the SPT data in this combination 
changes this result to the higher value of 0.90lg'23- A 
summary of the constraints with different combinations 
of data is provided in Table Q. 

Future CMB observations by the Planck satellite, as 
well as ground-based observatories are expected to im- 
prove this constraint dramatically over the next five 
years, and thus confirm or rule out this prediction. 



III. PRECISION TESTS OF GRAVITY 

Gravity on millimeter to solar system scales is well 
described by General Relativity, which has passed many 
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FIG. 2: Constraints at the 95% confidence level for Gn/Gr from WMAP 7-year (background, green), 
WMAP+ACT+SPT (middle, blue) and WMAP+ACT+SPT+Sne+BAO+Hubble data (front, red). The white lines 
show the 68% confidence levels. Note that the Gravitational Aether prediction is Gn/Gr = 0.75, while in General 

Relativity Gr = Gn ■ 




0.12 0.14 0.16 0.1S 0.20 0.22 0.24 0.26 0.6 0.7 O.S 0.9 1.0 1.1 0.1 0.2 0.3 0.4 0.5 



FIG. 3: Constraints at the 95% confidence level for Gn/Gr from WMAP+ACT+Sne+Ly-a (background, green) and 
WMAP+ACT+SPT+Sne+Ly-a (front, red). The white lines show the 68% confidence levels. Note that the 
Gravitational Aether prediction is Gn/Gr = 0.75, while in General Relativity Gr = Gn- 



precision tests on these scales with flying colors (see e.g., 
p3| for a review). That is why it is hard to imagine 
how an order unity change in the theory such as that 
of Eq. (0) can be consistent with these tests, without 
introducing any fine-tuned parameter. In this section, 
we argue that nearly all these tests are with gravitational 
sources that have negligible pressure or vorticity, which 
source deviations from GR predictions in gravitational 
aether theory. 



A. Parametrized Post-Newtonian (PPN) formalism 

In Sec. 0, we argued that for any perfect fluid with 
constant equation of state, w, the solutions of gravita- 
tional aether theory are identical to those of GR with 
a renormalized gravitational constant oc (1 + w). How- 
ever, for generic astrophysical applications, w is not con- 
stant except for pure radiation, or in the pressureless 
limit of w — 0. Focusing on the latter case, and given 
that pressure is 1st order in post-Newtonian expansion, 
we can quantify the gravitational aether theory through 
the Parametrized Post-Newtonian (PPN) formalism. 

The Parametrized Post-Newtonian (PPN) formalism 
is defined in a weak field, slow motion limit, and de- 
scribes the next-to-Newtonian order gravitational effects 
in terms of a standardized set of potentials and ten pa- 
rameters. These PPN parameters will be determined by 



solving the field equations ([l]) order-by-order with a per- 
fect fluid source in a standard coordinate gauge. The 
conventional introductory details of the formalism will 
be skipped over (see (34|] for a more detailed explanation 
of the procedure and the general PPN formalism). 

To be clear, though, we will assume a nearly globally 
Minkowskian coordinate system and basis with respect to 
which, at zeroth order, the metric is the Minkowski met- 
ric [Silv = Tjfiv) ancl the fluid four- velocity is purely 
timelike (u° — l,u l — 0). The stress-energy tensor is 
taken to have the form T^ v = (p + pll + p)u^u v + pg^ v 
where u M , p, II and p are the the unit four- velocity, 
rest-mass-energy density, internal energy density, and 
isotropic pressure of the fluid source, respectively. The 
fluid variables are assigned orders ofp^II^2^ u 2^ 

1PN. In the weak field limit, the metric can be written as 
a purturbation of the Minkowski metric: g^v = r\^ v 

V ' 

The components of the metric perturbations h^ iv with 
respect to this basis will be assumed to be of orders: 
h 0Q ~ 1PN + 2PN, hij ~ 1PN, and h oi ~ 1.5PN. This 
choice preserve the Newtonian limit while allowing one to 
determine just the first post- Newtonian corrections. Fur- 
thermore, the aether four-velocity u'^ will be assumed to 
be of the same order as that of the matter fluid. 

Solving (§) to 1PN gives p' = — p/4, which can be used 
in (0) to solve for 500 and gij to 1PN: 

h o - 2(7 (8) 
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2U6i 



(9) 



where U is the Newtonian potential and the follow- 
ing gauge condition is imposed: djhij — \(dihjj — 
dihoo)- Comparing the continuity equations for matter 
and aether (i.e. the "time" component of (||) to 1.5 PN), 
it can be shown that 



u 11 - u l = f 



(10) 



where t % satisfies V l (tip) = 0. This implies that the 
rotational component of aether is not fixed by matter 
within the PN expansion formalism. Here we will make 
the assumption that t % = so that aether is completely 
dragged by m atter. We will discuss this choice further in 
Section ( pIB| ). 

Previously we mentioned that in this case, an exact 
solution for u'^ and p' exists when matter has a constant 
equation of state. (It is worth noting that in the t % = 
case, higher PN equations appear to imply a nonstandard 
condition on the pressure V a (u a p) = 0, which is satisfied 
for a constant equation of state.) Using this solution and 
an additional gauge condition dih^i = SOqII, the field 
equations can be solved for goi and <?oo to 1.5PN and 
2PN, respectively: 



h 0i = - 7 -Vi - ~W U 



(11) 
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h 00 = 2U- 2U 2 + 40! + 40 2 + 203 + 6(1 + -)0 4 , (12) 

where Appendix ^ includes the definition for all poten- 
tials. Collecting all the results (§), (|), (|ll|), and @ in- 
dicates that all metric components are as in standard GR, 
except for the term in g 00 with the pressure-dependent 
potential £4- Consulting the parametrization rubric indi- 
cates that all PPN parameters have the standard values 
except £4, which equals 



<4 



1 

3' 



(13) 



which was already pointed out in |T(|. Notice that £4, i.e. 
the anomalous coupling of gravity to pressure is the only 
PPN parameter that is not measured experimentally, as 
one needs to probe the relationship between gravity and 
pressure of an object with near-relativistic pressures. A 
notable exception is observation of neutron stars (or their 
mergers, via gravitational wave observations), which can 
potentially measure £4, assuming that the uncertainties 
in nuclear equation of state are under control |35| . 



B. Gravitomagnetic Effect 

In the previous section, we showed that rotation of 
aether is not fixed by matter in the non-relativistic 
regime. We further assumed that aether rotates with 



matter. Here we will argue that observational bounds on 
the gravitomagnetic effect provide a mild prefernce for 
this assumption. 

Space-time around a rotating object with a weak grav- 
itational field, like Earth, can be described in terms of 
a set of potentials. With appropriate definitions, these 
potentials satisfy equations analogous to Maxwell's equa- 
tions [[36j . The gravitomagnetic effect describes the drag- 
ging of spacetime around a rotating object and can be 
quantified by a gravitomagnetic field B defined as: 



B 



3r(r-S)-Sr 2 



2r 5 



S l = 2G' e 



t 0k 



x' j T°£d 3 x'. 



(14) 



(15) 



where r is the position vector measured from the center 
of the object, e* fc is the three-dimensional Levi-Civita 
tensor, and is the RHS of the field equations (Q) p6| . 
The gravitomagnetic field causes the precession of the 
orbital angular momentum of a free falling test particle. 
The angular velocity of this precession is pq] 
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(16) 



If aether is irrotational, T c °g 
of linearized theory and since G — o" G tv ; we have: 



T ok to within the accuracy 



^aether — ^ GR - 



(17) 



Gravity Probe B (GP-B) is an experiment that mea- 
sures the precession rate < f2 > of four gyroscopes or- 
biting the Earth. Recently, GP-B reported a frame- 
dragging drift rate of —37.2 ± 7.2 mas/yr, to be com- 
pared with the GR prediction of —39.2 mas/yr ('mas' 
is milliarc-second) |3j. Laser ranging to the LAGEOS 
and LAGEOS II satellites also provides a measurement of 
the frame-dragging effect. The total uncertainty in this 
case is still being debated; with optimistic estimates of 
10% — 15% (e.g., fl38|)) and more conservative estimates 
as large as 20% - 30% (e.g., |§). 

Therefore, we conclude that even though perfect co- 
rotation of aether by matter is preferred by current tests 
of intrinsic gravitomagnetic effect, an irrotational aether 
is still consistent with present constraints at 2a level. 



IV. CONCLUSIONS AND DISCUSSIONS 

In the current work, we studied the phenomcnological 
implications of the gravitational aether theory, a modi- 
fication of Einstein's gravity that solves the old cosmo- 
logical constant problem at a semi-classical level. Wc 
showed that the deviations from General Relativity can 
only be significant in situations with relativistic pres- 
sure, or (potentially) relativistic vorticity. The most 
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prominent prediction of this theory is then that gravity 
should be 33% stronger in the cosmological radiation era 
than GR predictions. We showed that many cosmologi- 
cal observations, including CMB (with the exception of 
SPT), Ly-a forest, and 7 Li primordial abundance prefer 
this prediction, while deuterium may prefer GR values. 
We then examined the implications for precision tests of 
gravity using the PPN formalism, and showed that the 
only PPN parameter that deviates from its GR value is 
(4, the anomalous coupling to pressure, that has never 
been tested experimentally. Finally, we argued that cur- 
rent tests of Earth's gravitomagnetic effect mildly prefer 
a co-rotation of aether with matter, although they are 
consistent with an irrotational aether at 2a level. 

In our opinion, the fact that gravitational aether has 
the same number of free parameters as GR, and is yet 
(to our knowledge) consistent with all cosmological and 
precision tests of gravity at 2a level, indicates that this 
theory could be a strong contender for Einstein's theory 
of gravity. 

Future observations are expected to sharpen these dis- 
tinctions. In particular, the most clear test will come 
from the Planck CMB anisotropy power spectrum that 
is expected to be released in early 2013. Judging by 
the predictions for constraints on the effective number 
of neutrinos, Planck should be able to distinguish GR 
and Aether at close to 10er level |23| . 

Another interesting implication of this theory is for the 
cosmic baryon fraction. As we increase the gravity due 
to radiation (or effective number of neutrinos), we need 
to increase the dark matter density to keep the redshift 
of equality constant, since it is well constrained by CMB 
power spectrum (see e.g., [^l]). This implies that the 
total matter density should be bigger by a factor of 4/3 
(Fig. ||). Given that baryon density is insensitive to 
this change, the cosmic baryon fraction will decrease by 
a factor of 3/4, i.e. from 17% fl| to 13%. This could 
potentially resolve the missing baryon problem in galaxy 
clusters E(J, as well as the deficit in observed Sunyaev- 
Zel'dovich power spectra, in comparison with theoretical 
predictions j|| ||. 
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Appendix A: Aether perturbations through equality 

Here we present a consistent treatment of cosmologi- 
cal scalar perturbation theory for Gravitational Aether 
(GA). As we argued in Section y, when matter is ap- 
proximated by a perfect fluid with density p, pressure 
v = wp (w constant), and four velocity = f x , . 

r 1 v 11 V —ds 2 

(i.e. T M „ = (1 + w)pu^u lJ + wpg^y), u'^ = and 
p' = (i+^Xa™-!^ solve (|) and (|) and the GA field 
equation ([j]) becomes 

(Stt^G^ = G N (1 + w)T^ u . (Al) 

In cosmology, therefore, if the constituents of the uni- 
verse are matter and radiation and they are separately 
conserved, the GA field equations become 

4 

(8tt)~ G m „ = G N T™ V + -G N T^, (A2) 

where m and r stand for matter and radiation respec- 
tively. This approximation, of course, is false when in- 
homogeneities are considered since baryons and photons 
interact strongly. Therefore, we shall perturb about this 
exact solution and do a consistent treatment of cosmo- 
logical scalar perturbation theory. 

In what follows, b, dm, m, and r stand for baryon, dark 
matter, matter, and radiation respectively, and all barred 
quantities are unperturbed. Following p2], we will use 
the Conformal Newtonian Gauge: 

ds 2 = a 2 (r){-[l + 2TP{r,x)]dT 2 + [1 - 2<p(T,x)]dx i dx i }. 

(A3) 

To linear order in perturbation theory, the matter energy- 
momentum tensor takes the form 

T° = -ip + Sp) (A4) 

A = (P + P)— (A5) 
a 

T) = (p + 6p)5 ij +E i j , (A6) 

where Y?j is the traceless anisotropic shear stress pertur- 
bation and 

S P = P- P;5p = p l -p;8u^ = u^-u tJi , (A7) 

where i = {dm,b,r}. In our coordinate system u° = -, 
uq = —a, and Ui — u l — 0. The Aether pressure and 
four-velocity perturbations are defined as follows: 

p' = ~+ S P ', < = < m + Su,. (A8) 

Dark matter only interacts gravitationally and is sepa- 
rately conserved. We assume that there is negligible en- 
ergy transfer between baryons and relativistic particles 
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(i.e. V^fpftU^) = 0). Then, to first order in perturbation Conservation of G M „ and T^ v implies 
theory (Q) and (||) give: 



34 5p> = ^dtfu' 
a z 4 



Pb 



5w l ) 



di6p' 



4 a 



(A9) 
(A10) 



where 8w l 



5u l b = a 



l (8uf m - 5u\) and 5u l = 



Taking the comoving divergence of (A.10) and 



applying the comoving Laplacian to ( |A9|) , we can elimi- 
nate 5p' and get an equation for VL = diSu 1 : 



= (G R - G N )V V T. 



(A16) 



Defining p' 
tification u 



i(G R -G N ) 



and making the obvious iden- 



we see that this equation becomes 
exactly (||). Therefore, all of our solutions before can 
be used here after a trivial rescaling of the pressure. 
For example, if T^„ is a perfect fluid with equation of 
state w, exact solutions are obtained by u M = and 
P = (Gr — G;v)(l + w)(3w — l)p, which again just renor- 
malizes Newton's constant: 



3— d T (a 2 Vt) 



Pb 
ap m 



V 2 (6 b 



(All) 



where Sdn 



_ <V 



and we have used the 



this 



Pdm ' _ Pb 

fact that di8w l = ^-(5b — <Wi)- In Fourier space, 
equation can be numerically integrated for modes of dif- 
ferent wavelength, given the equations that govern 8 dm 
and 8b- Once f2 is known, (A9) can be used to find 
8p' . In the Conformal Newtonian Gauge, only scalar per- 
turbations are treated and we can ignore the rotational 
part of 8u z . This can also be physically motivated: let 
8u l — diUs + 08u l v where di6v,v = 0. Taking the curl 



of ( A10 ), it follows that V x Suy oc \. As a result, the 
rotational part of the Aether fluid decays and it doesn't 
play a major role in cosmology. As a result, given 57 we 
can find 8u l in Fourier space (dj — > ikj): 



k z 



G 



Gn — > G e ff(w) = Gn{3w—P- + (1 - 3iu)}. (A17) 

Note that G ef f(w = 0) = G N and G eff (w = 1/3) = Gr. 
Again, if matter and radiation are separately conserved 
in a cosmological setting, (A15) becomes 



- 1 s~i rpim 1 f > T^r 



(A18) 



More importantly, when Gr — Gn, these field equations 
reduce to those of General Relativity (GR) (this is 
true in the cosmological case because V^u^ ^ 0, which 
means that the conservation of Aether implies that its 
pressure vanishes identically). Also when Gr = |Gat, 
the GAA field equations reduce to those of GA, with the 
appropriate rescaling T' — j^—T^ u . Therefore, fitting 
this theory to data, we will be to make a likelihood 
plot of and see how far away the best fit is from the 



( Al2 ) GA and GR predictions. 



where k 2 = Sijkikj. Similarly, 

M = Him - H = % %( 8 ^ ~ (A13) 

To first order in perturbation theory, the GA field equa- 
tions now take the form 

(SttGjv)-^ = T™ + \t; v + (A14) 



with e 00 = 0, e 0i = ^-(Su, + j^Swi), and 



e ij = \a 2 8p'gij 



Because of the similarity of the underlying equa- 
tions, the linear perturbation theory of the GGA field 
equations is very close to those of GA, which we already 
described. We treat all matter perturbations as before 
and perturb similarly: 

p=(G N - G R ) Pm + Sp, = u-fT + Su^. ( A1 9) 
The equations of interest are (in Fourier space): 

3Hd T (a 2 n) + a(r)k 2 n = k 2 ^(8 dm - S b ) (A20) 

pmo 



Having both the left and right hand sides of this 
equation, we can now solve for the scalar perturbations. 
However, this does not provide an obvious way of 
checking the prediction of this theory, namely = |. 
This can be easily accommodated for by having field 
equations that contain Gr as a constant, and reduce 
to General Relativity and GA for Gr — Gn and 
Gr, = I Gat respectively. Consider the field equations 
(which we will refer to as the Generalized Gravitational 
Aether (GGA) field equations) 



in) 1 G llv \g tlv \ = GrT^ v + (Gn — GR)T^g^ u + T f 



T„ 



(A15) 



SP = ( G *- g">*" [ + J*L_(f 6 - S dm )] (A21) 
6tL apm 



k 



(A22) 



and we have recognized that 4 

6 p m fan 



where H = 

is fixed by the values at the present time. O nce ( A2C) 
i s sol ved for J7, Sp and 8u 3 are determined by (A.21) and 
( |A22| ), respectively. At long last, to linear order in per- 
turbation theory, the GAA field equations read 



(8tt) 1 G fliy — GnT™ + GrT^ u + e 



pu 



(A23) 
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where 
eoo = 



(A24) 



e«y = i^(G N - G R )(a 3 p m )[n + - 5, 



("Pino 



dm ) 



(A25) 



(Gr - Gn) " [0 H — r^—{h ~ Sd m ))Si 
3H ap mo 



(A26) 



Having both the left and right hand sides of ( A23 ), the 
scalar perturbations can now be consistently solved for. 



The potentials are all of the form 



F(x) = G N / d 6 y 



13.. Pfa)f 



\x - y\ 



and the correspondences F : / are given by 



(B4) 



Appendix B: PPN notations 

The metric components are in terms of particular po- 
tential functions, thus defining the PPN parameters: 



U 

4>w 



.goo = -l + 2J7-2 / 3[/ 2 -2^ + (27 + 2 + a 3 + Ci-2O0i A 
+ 2(3 7 - 2/3 + 1 + C 2 + 002 + 2(1 + Cs)03 
+ 2(3 7 + 3C 4 - 20^4 - (Ci ~ 2£)A (Bl) V t 

9ij = (l + 2 7 C/)% (B2) 

got = --(47 + 3 + ax -a 2 + Ci 
- 1(1 + 02 -Ci + 20Wi 



(B3) 



n 



04 : p/p 



(a - y)j r (y - z )j <f - z )j 



d 3 zp(z)- 

F - y| 2 L F _ z \ 

(vi(x - 7/),) 2 
\x - y\ 2 



W - z \ 



w, 



\x - y\ 2 



(B5) 
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